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We propose a simple Hamiltonian to describe the motion and the merging of Dirac points in the 
electronic spectrum of two-dimensional electrons. This merging is a topological transition which 
separates a semi-metallic phase with two Dirac cones from an insulating phase with a gap. We 
\ calculate the density of states and the specific heat. The spectrum in a magnetic field B is related 

to the resolution of a Schrodinger equation in a double well potential. They obey the general 
scaling law e„ oc B 2 ^ 3 f„(A/ B 2 ^ 3 . They evolve continuously from a sfnB to a linear (n + 1/2)B 
dependence, with a [{n + 1/2)B] 2 ^ 3 dependence at the transition. The spectrum in the vicinity of 
the topological transition is very well described by a semiclassical quantization rule. This model 
describes continuously the coupling between valleys associated with the two Dirac points, when 
approaching the transition. It is applied to the tight-binding model of graphene and its generalization 
£NJ . when one hopping parameter is varied. It remarkably reproduces the low field part of the Rammal- 

Hofstadter spectrum for the honeycomb lattice. 
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I. INTRODUCTION 



The main interest of graphene from the fondamental point of view is that the low energy electronic spectrum 
(around the band center e = 0) is linear, exhibiting the so-called Dirac spectrum around two special points K and 
K' at the corner of the Brillouin zone (BZ). 1 As a consequence, the density of states varies linearly with energy. In a 
magnetic field B, the energy levels around e = vary as e n (B) cx \J nB, with a two- fold degeneracy corresponding to 
the two valleys near K and K' . The considerable development of research on graphene is partly due to this unusual 
1 , ' spectrum. 2 

The electronic spectrum of graphene is very well described by a tight binding model on a honeycomb lattice, 
with three equal couplings t between nearest neighbors. 1 It has been soon realized that, by varying these hopping 
parameters, new interesting physics could emerge, in particular the existence of a topological transition separating a 
l— ~~ '■ metallic phase with two Dirac points and an insulating phase with a gap. 3-9 . When only one of the three hopping 
I \ parameters is modified (t' , see figure 10), the transition occurs when if = 2t. This model that we will call the t — t' 

■ model is generic and contains the essential physics of the more general case when the three hopping integrals are 
' different. 

Although such variation of hopping parameters may not be feasible in graphene, a transition could well be observed 
in other systems like the organic conductor a — (BEDT — TT F) 2 Iz W12 or an artificial lattice of cold atoms, 13-16 
, where the motion of Dirac points may be induced by changing the intensity of the laser fields. 
[**«- ■ The spectru m in a magnetic field of the honeycomb lattice has been first considered at low field by McClure 17 
who found a V nB dependence of the energy levels near the band center, and R. Rammal calculated the so-called 
" Hofstadter" spectrum which describes the fractal broadening of the Landau levels, due the the competition between 
magnetic field and lattice effects. 18,19 More recently the Hofstadter-Rammal spectrum has been studied when hopping 
parameters are modified, 20 and it was found that at the topological transition the low field spectrum exhibits a new 
dependence with the magnetic field of the form [(rt + l/2)£>] 2 / 3 . 4 This is due to the peculiar character of the dispersion 
relation at the transition: it is linear in one direction and quadratic in the other one. Quite recently it has been 

■ proposed that such a hybrid spectrum and the subsequent structure of the Landau levels may exist in VOijTiOi 
nanostructures. 21 

In a recent paper, we have studied under which general conditions a pair of Dirac points in the electronic spectrum 
of a two-dimensional crystal merges into a single point at the topological transition. 7 We have derived a low energy 
Hamiltonian that describes the physical properties near the transition. The present paper considers in details these 
physical properties. The main interest of this work is to present the simplest model which continuously describes 
the merging of two Dirac points and the subsequent gap opening, and to calculate several physical properties in the 
vicinity of the transition. We wish to stress that the interest of this work is not only to tune continuously the coupling 
between two Dirac valleys, but also to study the combination between two distinct and quite interesting dispersion 
relations : the linear dispersion relation and a dispersion relation with a saddle point. Indeed, the two Dirac points 
are always separated by a saddle point, and the merging of Dirac points is obviously accompanied by the merging with 
the saddle point. At this merging, a linear density of states characteristic of the Dirac point approaches a logarithmic 
density of states characteristic of a saddle point (in 2D). 



2 



The paper is organized as follows. In the next section, we consider a general tight binding problem in 2D, with 
two atoms per unit cell and study under which general conditions Dirac points may merge. In section III, we study 
several properties of this Hamiltonian. In particular, we relate the Landau level spectrum to a one-dimensional double 
well problem, where the two wells correspond to the two valleys around the Dirac points. When approaching the 
transition, the potential well vanishes and the potential becomes quartic at the transition. In section IV, we show 
how to relate the parameters of the universal Hamiltonian to the parameters t and t' of the tight binding model on 
the honeycomb lattice, the so-called t — t' model. In this way, we reproduce with a very good accuracy the low field 
part of the butterfly spectrum spectrum and the lifting of the valley degeneracy of the Landau levels 



II. CONSTRUCTION OF THE UNIVERSAL HAMILTONIAN 



We consider a two-band Hamiltonian for a 2D crystal with two atoms A and B per unit cell. This two-band 
Hamiltonian is naturally given in terms of the 2x2 matrix 

H(k) = ( h h AA fJ h h AB % ) , 
\ hBA(k) h BB {k) J 

with the 2D wave vector k. Time-reversal symmetry (H(k) = H*(—k)) imposes fiAB(k) = h,BA*(k) = f(k) and, 
together with hermiticity, real symmetric diagonal terms h,AA(k) = h-AA(-k) (h B B(k) = /ibb(— k)). Furthermore, we 
consider a 2D lattice with inversion symmetry such that JiAA(k) = h B B(k). The resulting energy dispersion reads 
e±(k) = hAA(k) ± |/(fc)|, and we will even set IiAA(k) = because this term simply shifts the energy as a function of 
the wave vector but does not affect the topological properties of the semi-metal-insulator phase transition discussed 
below. 

We, therefore, discuss from now on the Hamiltonian in its reduced form 



where the off-diagonal terms have the periodicity of the Bravais lattice and may be written quite generally in the 
form: 



/(fc) = X> m „e- i£ -*™» , (2) 

m,n 

where the t mn 's are real, a consequence of time-reversal symmetry 7i(k) = H*(—k), and R mn = may +na 2 are vectors 
of the underlying Bravais lattice. 

The energy spectrum is given by e(k) — ±|/(fc)|, and the Dirac points, that we name D and — D are solutions of 
f(D) = 0. Since f(k) = f*(—k), the Dirac points, when they exist, necessarily come in by pairs. 22 The position D of 
the Dirac points can be anywhere in the BZ and move upon variation of the band parameters t mn . Around the Dirac 
points ±L>, the function f(k) varies linearly. Writing k — ±D + q, we find 

f(±D + q) = q-(±v 1 -lv 2 ) (3) 

where the velocities v\ and vi are given by 

V\ — ^ ^ tmnBmn sin D • Rmn 
mn 

V2 = ^t mn Rmn COS D ■ R mn (4) 

mn 

Upon variation of the band parameters, the two Dirac points may approach each other and merge into a single 
point Do. This happens when D = —D modulo a reciprocal lattice vector G = pal + qa'2, where a\ and a% span the 
reciprocal lattice. Therefore, the location of this merging point is simply Dq = G/2. There are then four possible 
inequivalent points the coordinates of which are D n = (pal + qa%)/2, with (p, q) = (0, 0), (1, 0), (0, 1), and (1, 1). The 
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condition /(-Do) — Emnf"^'""'™ = 0, where j3 mn = V m + 9 n ! defines a manifold in the space of band parameters. 
As we discuss below, this manifold separates a semi-metallic phase with two Dirac cones and a band insulator. 

In the vicinity of the D point, / is purely imaginary (vf = 0), since sin(G • R mn /2) = 0. Consequently, to lowest 
order, the linearized Hamiltonian reduces to TH. = q • v%<r v ', where v% = ^2 mn {—l)^ mn tm.nRm.n- We choose the local 



reference system such that v% = c y y defines the y-direction. 23 In order to account for the dispersion in the local 
^-direction, we have to expand f(Do + q) to second order in q: 

f(D + q) = -iq- £(-l)^»i mn (g • R mn f . 

mn 

Keeping the quadratic term in q xi the new Hamiltonian may be written as 



(5) 



H (q) 

where the effective mass m* is defined by 





2m* lc y1y 



2m* 



ic y q y 
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£(-i)' 



u mn r ^ mn ,x ' 



(6) 



(7) 



and where R m n,x is the component of R mn along the local x-axis (perpendicular to w>). The terms of order q y and 
q x q y are neglected at low energy. The diagonalization of Ho(q) is straightforward and the energy spectrum 



e = ± 



2 2 , 
C y q y + 



2m" 



1/2 



(8) 



has a remarkable structure: it is linear in one direction and quadratic in the other. From the linear-quadratic spectrum 
which defines a velocity c y and a mass m* , one may identify a characteristic energy : 



Emn( jj m "t mn R mn ] 2 
Emni-^-^tmnR* 



(9) 



The merging of the Dirac points in Do marks the transition between a semi-metallic phase and an insulating phase. 
In this paper, we concentrate on the properties of the spectrum in the vicinity of the merging. The transition is driven 
by the parameter 



A = /(Do) = J2(-l) 0mn t„ 



(10) 



which changes its sign at the transition. This parameter A therefore drives the transition. In the vicinity of the 
transition, the Hamiltonian has the form 



H(q) 



with the spectrum 



2m* 



"yly 



A + 2^ + ic yly 



±\/(A + 



2m 



■f + qlcl 



(11) 



(12) 



The Hamiltonian (11) has a remarkable structure and describes properly the vicinity of the topological transition, 
as shown on Fig. 1. When m*A is negative (we choose to* > without loss of generality), the spectrum exhibits the 
two Dirac cones and a saddle point in D (at half distance between the two Dirac points). Increasing A from negative 
to positive values, the saddle point evolves into the hybrid point at the transition (A = 0) before a gap 2A > opens. 

In this paper, we study the spectral properties around the merging, in particular in the presence of a magnetic field. 
Moreover, we stress that this Hamiltonian has the general structure to describe the physics of Dirac points, even far 
from the transition, since it captures quite simply the coupling between the two valleys associated with the two Dirac 
points. In particular, we can relate the coupling between valleys to a double well potential problem. For this reason 
we name it a universal Hamiltonian. 
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qo, m* 


c x = qD/m* 


A = 


-q 2 D /2m* 


qD, c x 


m* = q D /c x 


A = 


-c x q D /2 


m*, c x 


q D = m*c x 


A = 


-m*c 2 /2 


m*, A 


c x = v / ~2A/m* 


qo = 


V-2m*A 


Cx, A 


qo = -2A/c x 


m* = 


= -2A/c 2 


qo, A 


m* = -q%/2A 


c x = 


-2A/q D 



TABLE I: In addition to the velocity c y , the universal Hamiltonian is described by two independent parameters (left column) 
from which two other parameters may be deduced 



III. PROPERTIES OF THE UNIVERSAL HAMILTONIAN 




FIG. 1: Evolution of the spectrum when the quantity A is varied and changes in sign at the topological transition (arbitrary 
units). The low-energy spectrum stays linear in the q y direction. 

Without loss of generality, we assume to* > 0. When A varies from negative to positive values, a topological 
transition from a semi-metallic phase with two Dirac cones and a band insulator with a gapped spectrum occurs. At 
the transition, the spectrum is hybrid, a reminiscence of the saddle point in the semi-metallic phase, see figure (1). 

When A < 0, the spectrum exhibits two Dirac points the position of which along the x axis is given by iqD with 



q D = V-2m*A (13) 
and the linear spectrum around these Dirac points is characterized by the velocity c x along the x direction : 

qD (M\ 

to* V m* 

The two Dirac points are separated by a saddle point at position qs — whose energy is ±|A|. The mass to* describes 
the curvature of the spectrum at the vicinity of this point along the x direction. When A < varies, the Dirac points 
move along the q x axis. 

The energy dispersion relation (12) is characterized by three parameters, the velocity c y along the q y direction, the 
mass m* along the q x direction and the gap A. Alternatively, it can be characterized by the distance 2qu between 
the Dirac cones or the velocity c x , or by any combination of two among the four parameters to*, A, c x , or qo- In table 
(I), we write explicitly all the combinations between these parameters. The universal Hamiltonian describes properly 
the vicinity of the Dirac points. The spectrum can be linearized along the q x direction, to recover a anisotropic Dirac 
equation in the vicinity of each point with a velocity c x given by (14) 
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n = I ±c x 5q x - ic y q y \ 

\ ±c x 8q x + ic y q y J 

where Sq x = q x - q D . 

In section IV, we discuss which combination of parameters should be taken in order to properly describe the low 
energy physics of the t — t' model of the honeycomb lattice. 



A. Density of states 

We have calculated the density of states for the energy dispersion (12). We find 

where K{x) is the complete elliptic integral of the first kind. 24 In the low energy limit, one recovers the familiar linear 
energy dependence p(e) = 1 - e. The density of states exhibits a logarithmic divergence at e = |A|, due to the saddle 
point. It is plotted in figure (2) for a fixed mass to*, and upon variation of the parameter A. When approaching the 
transition, the weight of the logarithmic singularity vanishes and, at the transition, one recovers the density of states 
found in ref. 4, given by 



Cy 



(16) 



where C = ^K{\/y/2) = r(l/4) 2 /(47r 5 / 2 ) ~ 0.188. 

Above the transition, there is a finite gap A > and the density of states has a jump at e = A 



e < A 
e > A 



p(e) = 
P(e) = 




(17) 



B. Specific heat 



Since the spectrum has the electron- hole symmetry p(e) = p(—e), the chemical potential is temperature independent 
and fixed at p = for the undoped system. Therefore the specific heat has the general form 



where (3 = l/(fcsT). Using the expression (15) of the density of states, we obtain, on the metallic side (S < 0) 



where the function f(T/\A\) is given by 



C (T)=4k B ^(k B T)^f[^-) (19) 
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FIG. 2: Evolution of the density of states on the metallic side of the transition (A < 0), plotted for A = —1 and A = —0.2. 
For finite A, there is a linear dependence at low energy <C |A|, followed by a logarithmic singularity at |A|. At the transition 
the density of states varies as y/e (dashed curve). Above the transition, there is a finite gap (A > 0) and the density of states 
has a discontinuity (dotted curve for A — 0.7). 



1/24 2V~tx 3 „, I Ax s , r * 5/2 „, I2xt+1 



/(*)=/ , o~KU- -)dx+ K(\ — )dx (20) 

It is plotted on figure (3) and it has the following limits 
t -> /(*) -» ^C(3)Vt ^ 4.248Vi 

t - oo /(t) - ^(V8- l)r(l/4) 2 C(5/2) ~ 1.889 

so that the specific heat interpolates from a T 2 behavior far from the transition to a T 3 / 2 behavior at the transition. 
Similarly, above the transition, in the insulating phase (S > 0), we find, using (17): 



tt c v 

where the function g(T/A) is given by 



C(T) = 4k B ^(k B T)^g(l) (21) 



x 



5 / 2 l2xt- 1 



9(f) = —^- K (\— A ) dx ( 22 ) 

Ji/2t cosh x V 4a; 

The function g(t) varies exponentially at small t, so that the specific interpolates from the T 3 / 2 behavior at the 
transition to an activated behavior. 



C. Landau quantization 

We now calculate the spectrum of the universal Hamiltonian in the presence of a magnetic field B. In the Landau 
gauge A = (0, Bx, 0), the substitution q x —fq x — eBy leads to the new Hamiltonian 

H= ( A + ±TO*w 2 y 2 - ic y q y \ 

y A + ±m*cj 2 y 2 - ic y q y J 

where, as usual, y = y — q x /eB. The spectrum in the bulk does not depend on the shift q x /eB. We introduce the 
dimensionless variables Y = y/a and q y = /aP, so that [X, P] = i. We obtain the effective magnetic length for 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 



FIG. 3: Top figure: function fit — T/\A\) entering the expression of the specific heat (19), in the metallic phase (A < 0). 
Bottom figure: function g(t = T/A) entering the expression of the specific heat (21), in the insulating phase (A > 0). 



1/3 

this problem a = ^ m*L ,A ) = (^ m * c y^B) 1 ^ 3 ^ where £ = \j\feB is the usual magnetic length. The Hamiltonian is 



rewritten as 



n 



m*uj 2 c c 2 y 



1/3 



S + Y 2 -iP 

5 + Y 2 + iP 



(23) 



Squaring this Hamiltonian, we have to solve the effective Schrodinger equation 



2/3 



* 2 2 

m W c °y \ ( 2 , ts , v-2\2 „-„rD t^2i\ ,,.A,B 



(P 2 + {5 + Y 2 ) 2 -is[P,Y 2 ]) ^ 



where s = ±1 corresponds to the two sites A and B. We have introduced the dimensionless gap 



We now have to diagonalize the effective hamiltonian 7i e // 



2 C 2 \ 1/3 B 2/3 



(24) 



Heff = P +{S + Y 2 Y ~2sY 



(25) 



and the eigenvalues e„ of the original problem (23) are related to the eigenvalues E n of this effective Hamiltonian (25) 
by 



± 



■m*u 2 c c 2 v 



1/3 



A 



^E n {8) = ±j^E n {8) 



(26) 



where 8 is given by (24). We thus obtain the general scaling behavior of the Landau levels. 




e n oc B 2 / 3 /n(A/B 2/3 ) (27) 

When varying 5, this Hamiltonian has the remarkable property to describe continuously the Landau level spectrum 
from the e„ oc \fnB dependence with double degeneracy for well separated Dirac cones to the e„ oc (n + 1/2)5 usual 
dependence for a massive particle. The physics behind is that for negative 5, the problem is similar to the one of 
a particle in a double well potential. In the limit of large negative 5, that is far from the transition or in a weak 
magnetic field, the potential has two well separated valleys which are almost uncoupled. This corresponds to the 
situation of two independent valleys. Note that in this limit the energy shift between the two valleys is 2\/~8. When 5 
diminishes, we progressively increase the coupling between valleys. The degeneracy of Landau levels is progressively 
lifted, as shown on Fig. 4. 

We have solved numerically the Hamiltonian Ti e ff- Eigenvalues are given on Fig. (5) as functions of 5. We now 
comment our results and the different limits. 

• If S < 0, we have to solve the problem of a double well potential (Figure 4. a) (and independent wells in the limit 
\S\ ^> 1). The potential has two minima for Yq = ±y/\5\. In real space, the distance 2Yq corresponds to the distance 
2y = 2aY Q = 2q D £ 2 B , where t B is the magnetic length. An expansion around these minima x gives the 

effective Hamiltonian : 

H e// =P 2 + 4|% 2 T 2 S /H . 
Introducing the new variables p' = P/(v / 2|<5| 1 / 4 ) and x' — xy^lS] 1 ^ 4 , this Hamiltonian reduces to 



H eff = 2^W\(p>*+ x > 2 ±l) 
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-6 -4 -2 2 4 



FIG. 5: Energy levels E n {5) of the Hamiltonian Ti e ff — P 2 + (5 + Y 2 ) 2 — 2Y . We have plotted asymptotic analytical behaviors, 
E„ — Any/— 5 for large negative S and E n — 5 2 + 4(n + 1/2) Vo for large positive 5. . 




-4 -2 2 4 



FIG. 6: Energy levels e' n (5) = e n (S)/(m*u> 2 Cy/2) 1 ^ 3 as a function of the dimensionless gap S. The dots on the 5 = axis 
indicate the semiclassical levels of the quartic Hamiltonian. 

with eigenvalues 

E„ = 4n^/\S\ (28) 
Each level n ^ is doubly degenerate, due to the twofold structure of the potential well. We deduce that 

a _ , A 2 

4 | (5 |3/2 n 

which can be written in the usual form, introducing c x = y/—2Ajm* 

e n = ±V2nec 2 B 

with the velocity c defined as 



C \J C X Cy 
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Each energy level e n is doubly degenerate. We recover the well-known result for two independent Dirac valleys, 
generalized here to the anisotropic case. 

• When 1 5 | diminishes, the potential barrier between the two valleys decreases and tunneling between the valleys 
removes the twofold degeneracy of each level (Figure 4.b). We can estimate the shift of the levels due to a finite A. 
The shift is proportional to the probability to tunnel between the two valleys. It scales as 5E n oc e -VV-E n d wnere 
the potential height V is proportional to \S\ 2 and the distance between valleys d is proportional to \/\6\- As a result, 
the level degeneracy is lifted as 



e -\s\^ „ e -#|A|VV* . (29 ) 

• At the transition point, 6 = 0, the energy levels are those of a modified quartic oscillator with a potential 
V(Y) = Y 4 — 2Y and they have been obtained in Ref. 4 and are well approximated by: : 

E n = C(n + 1/2) 4 / 3 (30) 
with C = 7r 2 [3v / 2/r(l/4) 2 ] 4 / 3 ~ 2.185. From eq. (26), we deduce the following dependence of the Landau levels 

e„ = ±A{m*cl) l l%n + l/2)c c ] 4 / 3 (31) 



with A = y/C/2 3 / 2 ~ 1.173. In ref. 4, we have studied in details the effect of the linear term in the potential V(Y), 
which only slightly change the above result. We have attributed the phase term 1/2 to the annihilation of the Berry 
phases attached to each Dirac point at their merging. This is also briefly discussed in the next subsection of this 
paper. 

• For large 8 > 0, the Hamiltonian can be expanded and transformed into a quadratic Hamiltonian 

Heft =P 2 + 5 2 + 25Y 2 -2Y~P 2 + S 2 + 26(Y - l/{25)) 2 (32) 
so that the spectrum is again the one of an harmonic oscillator 

E n = 5 2 + 2V2(n' + l/2)VS (33) 
and we recover a usual Landau spectrum in the gapped phase. 



£ „ = ±(A + ^ — (n+l/2K) (34) 
D. Berry's phase 

We now briefly turn to the structure of the wave functions, solutions of the universal Hamiltonian (11). They are 
of the form 




^=^( e ««r (35) 

where the two components refer to the two sublattices A and B. The phase 9^ is given by 

tang ? = CyQv 2 . (36) 



Note that the two valleys, centered on q x = ±qo = 2m* A, are described by the same wave function. The q 

dependence of 9^ is shown on Fig. 7 and exhibits a vortex structure around the two Dirac points. Each point is 
characterized by a Berry phase \§^9 q ■ dq = ±ir. Fig. 7 shows the annihilation of the two Berry phases at the 
topological transition. This is the reason why the Landau levels acquire a n + 1/2 dependence near and above the 
transition (see next section). 25 
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FIG. 7: Relative phase 9q of the two- component wave function. The parameter are chosen in arbitrary units m* = c y = 1. The 
four plots correspond respectively from left to right and then from top to bottom: A = — 1, — .3, 0, 1. In the insulating phase, two 
opposite Berry phases are attached to the two Dirac points. The Berry phases annihilate at the transition point. 




FIG. 8: Semiclassical quantization of area. When e < —A, the quantization of energy levels results from the quantization of 
orbits in each valley S(e)£% = 2imeB and the spectrum has the double valley degeneracy. When e > —A, above the saddle 
point, the quantization implies larger orbits which encircle the two Dirac valleys, and it reads S(e)£% = 2w(n' + l/2)eB . 



E. Semiclassical quantization and integrated density of states 

It is instructive to derive the energy levels from semiclassical Bohr-Sommerfeld quantization: along one period of 
the motion, the action must be quantized. This condition can be written as 

5(e) = 2^(n + 7 )^| , (37) 

where 5(e) is the area of a cyclotron orbit of energy e is reciprocal space. It is simply 5(e) = 47r 2 iV(e) where N(e) 
is the integrated density of states which can be obtained from expressions (15). The phase mismatch 7 is the sum 
of two contributions 7 = + 1b where 7a/ is the Maslov contribution and 7^ results from the Berry phase. We 
obtain, for e < —A: 

< e+a >*(\/S)- aE (\/S)] < 38 > 

where K(x) and E(x) are respectively complete elliptic integrals of the first and of the second kind. 24 This quantity 
represents the area enclosed by each of the two degenerate equal energy lines encircling one Dirac point (Fig. 8). 
The phase mismatch cancels here due to a finite Berry phase 73 = ±1/2, 25 so that the quantization condition is 
5(e) = 2imeB. 



5(e) = 



4\/2 y/m*{e- A) 
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Similarly, for e > — A: 



5(e) = 



(e + A)K ^ 



2c 



2AE 



2c 



(39) 



and the quantization condition is now S(c) = 2ir(n / + l/2)ei3. The contribution is canceled since the semiclassical 
trajectories enclose the two Dirac points and the Berry phase is 0. 

Figure (9) compares the real spectrum with the above semiclassical quantization. The approximation works very 
well except in the vicinity of the transition line e n = — A which corresponds to the energy of the saddle point. It is 
worth stressing that the semiclassical approximation describes perfectly well the vicinity of the topological transition 
(near 6 = axis on Fig. 9. The energy levels are given by the dimensionless equations 



with 



c < -A 
c > -A 



F_(r) = Vr - 1 
F+(r) = ^ 



{r + l)K 



e 
A 

(i 



I 2r 
r - 1 



37r n 

3tt n' + 1/2 
2^2 <5 3/2 



- E 




(40) 
(41) 

(42) 
(43) 



5 
4 
3 




-4 



-2 







2 



FIG. 9: Positive exact energy levels e' n (S) = c„(8)/(m*0JcCy/2) 1 ^'' i (dots) compared with the result of semiclassical quantization 
(lines). The dashed line c n = — S corresponds to e n = —A, that is to the energy of the saddle point. The discontinuity is due 
to the doubling of the area S(e) and to the cancelation of the Berry phase when energy crosses the saddle point. 



IV. APPLICATION TO GRAPHENE AND THE HONEYCOMB LATTICE 



We now propose that the effective Hamiltonian constitutes an excellent description of the low energy physics of the 
t — t' model including the two valleys. First we briefly recall the electronic structure of graphene, and assume more 
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generally that one of the three hopping parameters if between nearest neighbors may be different from the two others 
t, as shown on Fig. (10). 




FIG. 10: Honeycomb lattice with hopping integrals t and t' , and elementary vectors a\ and 3,2 discussed in the text. 



The tight-binding Hamiltonian couples sites of different sublattices named A and B. The eigenvectors are Bloch 
waves of the form 



\k) 



1 



(44) 



where \Rf), \Rf) are atomic states. The sum runs over vectors of the Bravais lattice. The Hamiltonian has the form 



(1), with 



f(k) =t' + te tk - ai +te 



(45) 



where a\ — a(-^, §), a-2 = a(~-^j |) are elementary vectors of the Bravais lattice, a is the interatomic distance, and 
t, if are shown in Fig. 10. In Cartesian units 



f(k) =t' + 2t cos yfc lB e l '^ a 



(46) 



The energy, given by e(k) — ±|/(fc)|, is shown in Figure (11) in the form of equal energy lines. 

The evolution of the low energy spectrum when if varies is plotted on Fig. 12. It is obviously well described by our 
universal Hamiltonian (compare with Fig. 1). In the following, we carefully map the two models on each other. 

When if — t, the energy vanishes at the two points D and D' located at the corners K and K' of the Brillouin zone 
K = 2a|/3 + 01/3, K' = af/3 + 2a2/3, where a\ and are reciprocal lattice vectors) or, in Cartesians units 



D = K 



( 2tt 2tt 



\3V3a 3a 



\3V3a' 3a / 

As if increases, the two points D and D' approach each other. Their position is given by 



/ 2 /4i 2 2tt 



(47) 



(48) 



They merge into the single point Do = (a| + a%)/2 = (0, |^) when t' = 2t (for t' > 2t, a gap opens between the two 
subbands). (a = 1 for shorter notations) 

We now concentrate on the vicinity of the KK' axis , that is the line k y = 2tt/3. An expansion near this line, gives 
(k y = 2-7I-/3 + q y ): 



f(k) = if - 2tcos^-fc 2 



3itq y cos 



V3, 



(49) 



-3-2-10 1 2 3 
FIG. 11: Isoenergy lines, fort'/t — 1, 1.5,2. 



t'=t 7 T y r=i.5f 





FIG. 12: Evolution of the low energy spectrum when t! approaches t. This evolution is very well described by the universal 
Hamiltonian (see figure (1). 



We now wish to describe this Hamiltonian by the universal Hamiltonian (11), that is 
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f(k)=A-ic y q y + ^ (50) 

for which we recall that fixing A and m* imposes the position ±qo of the Dirac points and the velocity c x (see table 
I). We are now facing several possible choices to properly introduce the effective Hamiltonian. We may choose to fix 
the mass m* and the parameter A by comparing the expansion of (49) near q x = 0: 



with (50). This leads to 



f(k) = t ' -2t+- A tq 2 x -mq y (51) 



A = t' - 2t , m* = ^ , c y = 3t (52) 



and qo and c x are obtained from table (I) and are plotted in Fig. 13. This is not a good choice because, if it properly 
describes the spectrum near q x = 0, it does not correctly describe the vicinity of the Dirac points iqD- 

We may also choose to fix A and the distance 2qo between the Dirac points and the velocity c y around the Dirac 
points 



2 4t 2 3 

A = t' - 2t , q D = —j= arctan y - 1 , c y = -t' (53) 

so that the mass and the velocity c x are deduced from table (I) and are plotted on Fig. 13. With this choice the 
velocity near the Dirac points is not correct, so that the low energy spectrum when the Dirac points are far apart 
cannot be reproduced. 

Among other possibilities we finally choose to fix A and the velocities c x and c y . Comparing (50) with the linear 
expansion of (49) near the Dirac points 

f(k) = ^it'q y ± v/3(* 2 - t' 2 /A)q x 
where the ± sign denotes the vicinity of the two points D and D' , we are led to choose the combination of parameters: 

of 

A = t'-2t , c x = y/3y/t 2 - 1' 2 /4 , c y = — (54) 
from which we deduce the effective mass 



m * = ^ = W^V) (55) 

This last choice of parameter is the best one since it properly describes the low energy spectrum with the correct 
velocities (see Fig. 13). Note that the low energy spectrum is not monotonic when if increases since the product 
c 2 = c x c y first increases and then decreases with t' (figure 14). 

Using these parameters, we can now apply the results of the universal Hamiltonian (figures 5, 6) to the specific 
t — t' model. We first introduce the reduced flux / = 4>/4>o, where (f> is the flux through one elementary cell of the 

honeycomb lattice, and (fro = h/e = 2n/e is the flux quantum. We have / = Ba 2 3\/3/(2(/) a ) — ^^-eB, since we have 
chosen a = 1, h = 1. From our study of the universal Hamiltonian, the energy levels are given by (26) 

t' - 2t 



e n (/) = y/E^S) (56) 

where the function E n (5) has been studied in section III. 3 (figure 5) and the parameter S is deduced from the 
parameters (54, 55). We have 
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[ m 




FIG. 13: Plots of the dependence of the quantities c x , m* and qo as a function oft' . The three lines correspond to three possible 
fits, where the mass m* is fixed (dashed lines), the position qD of the Dirac points is fixed (thin lines), or the velocity c x is 
fixed, as chosen in the text (thick lines). The dot indicates the correct variation in the t,t' model. 




FIG. 14: Non-monotonic behavior of the velocity c = ^/c x c y = ^—\ftJ (j±t? — t' 2 ) 1 ^ 4 as a function oft'/t. The velocity is plotted 
in units of ta/h. 



'2\ 2/3 t'-2t 1 
\n) [(2t + ^' 2 ] 1/3 / 2/3 ' 



In particular, in low field : 



e„ = ±yj2nec x c v B = ^2irnt'(4i 2 - t' 2 ) 1 / 2 ^/] . (58) 

Fig. 15 represents the energy levels for the honeycomb lattice with t' = 1.5t. The spectrum is represented as a 
function of the reduced flux /. In low field, the levels e have a y/nj behavior. Then the degeneracy of the levels 
is lifted as predicted in eq. (29), that is Ae cx e~#( 2t ~* ') 3/2 /^. The overall spectrum is quite well described by the 
semiclassical quantization rule explicited in section III. 5. 
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0.00 0.01 0.02 0.03 0.04 

f 

FIG. 15: Color on line. Red dots : energy levels of the tight binding model on the honeycomb lattice with t' = 1.5t. B0 Violet dots 
: energy levels calculated from the solutions of the ID Schrddinger equation with a double well potential V(Y) — (<5 — Y 2 ) 2 — 2Y . 
Full continuous lines : result of the semiclassical quantization rule. Dashed line : line e = —A = It — t . 



V. SUMMARY 



We have shown that the motion and merging of Dirac points in a two-dimensional crystal can be fully described by 
a simple 2x2 Hamiltonian with a linear dispersion relation in one direction, a massive term in the other direction, 
and gap term A. By varying A, a topological transition is driven, separating a semi- metallic phase with two Dirac 
points and a gapped phase. We have calculated analytically several quantities, such as the density of states, the 
specific heat, and the integrate density of states related to the area of semiclassical orbits. From this quantity, we 
obtain a simple semiclassical description of the Landau levels spectrum in a magnetic field B. More quantitatively, 
the problem in a magnetic field is related to a one-dimensional Schrodinger equation with a double well potential, 
whose potential barrier depends on the parameter A and the magnetic field B as A/B 3 / 2 . The spectrum of Landau 
levels scales as e n oc B 2 / 3 f n (A/ B 3 / 2 ). In the vicinity of the topological transition, it is very well described by the 
Bohr-Sommerfeld quantization rule S(e)l 2 B — 2n(n + 1/2). 

This Hamiltonian is appropriate to describe continuously the coupling between valleys which is usually neglected in 
the case of graphene, but which becomes important when approaching the topological transition. During completion 
of this paper, we have been aware of similar results in Ref. 26. 
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